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this note, we study the models proposed by Patel. We also study the large 



\^ ' Nc limits of Patel's model. To discuss the validity of the mean field theory 

results, the one-loop Coleman- Weinberg effective potential is calculated for 
Q^' Nc = 3. We argue that the quantum corrections vanish at large Nc when the 



energy of the so-called baryonic vertices scale with Nc- 



X 

c3 ■ I. INTRODUCTION 



The structure of the quantum chromodynamics (QCD) the vacuum has been a subject of 
intense investigations over the last two decades as qualitative understanding of hadron 
spectroscopy depends on the vacuum configuration. Since the QCD vacuum confines color 
degrees of freedom, it is highly non-trivial compared to the vacuum for QED. However, to 
address this issue in a sensible way from the fundamentals of QCD is a rather challenging 
task as the theory becomes strongly non-perturbative in the infrared region. 

It is customary to use effective Lagrangians in the low energy domain to describe the 
hadronic dynamics using the global symmetries of the theory. However, these models are 
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inadequate for studying the vacuum structure since the order parameters - the mesonic 
fields, are color singlets. 

Therefore, one has to take a path somewhere in between, explicitly keeping some of 
the colored degrees of freedom, like quarks, while truncating the other degrees of freedom, 
namely the gluons. These truncations are performed in various disguises - via potentials, 
bags, strings etc. One can also argue in favor of the roles played by topologically non- 
trivial classical solutions like instantons and subsequently perform a semiclassical analysis 
around these solutions. 

An economical picture for confinement is provided by the Meissner effect in type-II 
superconductors. It is argued that the behavior of QCD vacuum is dual to the type-II 
superconductors 0,^. In this picture, the chromoelectric field is screened in the "dual 
superconductor" just like magnetic field in ordinary superconductors. Therefore the chro- 
moelectric field between a quark and an anti-quark gets squeezed in a tube-shaped region 
like the Abrikosov fiux-tubes in type-II superconductors. This stringy picture for mesons 
dates back to the Nielsen-Olesen model [^], where the existence of a scalar field is crucial. 
But, apparently in QCD no such scalars exist. According to 't Hooft the role of the scalars 
can be mocked up in QCD via a suitable gauge choice and the so-called "Abelian Projec- 
tion" 0. Evidences for this type of string models come from the Regge behavior and the 
hadronization processes in high energy hadronic reactions where it gives a rather consistent 
picture 0. 

It is possible to incorporate baryons in this picture as well . To do so, one has to include 
special configurations called the "Y- vertices" , which allows A^^^ colored quarks to bind and 
hence form baryons. These vertices are exclusive to the models with Nc > 2 and the fiux- 
strings can form branches in these cases. Evidence for appearance of such objects in the 
magnetic condensation picture has also been argued by using Makeenko-Migdal equations 
p!0| . It has also been argued that these vertices are gluonic in nature and trace the baryon 



number [Oj. One is allowed to associate some energy with these vertices which is not 



associated with the fiux-tube itself. 



On the other hand, at high temperatures it has been argued from general principles 
|jl^,[13| that a deconfining phase in QCD exists. An important issue is then the nature of the 
phase transition. The order of phase transition can be found by using universality arguments 
and looking at appropriate spin models which belong to the same universality class as QCD 



for various color groups. In these spin models, the phase transition proceeds via the 
formation of domains at lower temperatures by cooling the system from a high temperature 
phase where one is guaranteed a disordered state which can be understood as the deconfined 
phase. 

One can use the flux-string models to study the same problem, though here one adopts 
a bottom-up approach - i.e. we start at lower temperature in the confining phase and then 
raise the temperature of the system. This forces the flux-string to fluctuate and at high 
temperatures the string will elongate and make a criss-crossed mesh. At the deconfining 
point this net will permeate the whole space. This way of addressing the deconfinement 



transition in a flux-string model was studied by Patel p5|J]i 

In the original work of Patel [|I^], a flux-string model was proposed on the lattice where 
the strings lived on the links and the quarks and the baryonic vertices were defined on the 
lattice sites subject to the Gauss law for the conservation of U{1) color charges modulo Nc- 
The model is simple enough to be studied analytically yet it retains much of the relevant 
physics in a model. Patel used mean field theory to evaluate the free energy of the system. 
In this treatment, the system showed phase-transition of second order for = 2 and of 
first order when A^^^ = 3 provided one takes into account the baryonic "Y- vertices" . This 
result qualitatively agrees with the result found by looking at the center of the gauge group 



1^ . The difference between Nc = 2 and A^c = 3 arises from the presence of the so-called 
"Y-vertices" in the string picture. In fact, this is consistent with the model of Patel. 

In fact, the gluonic vertices resolve the following puzzle. In absence of these vertices and 
quarks, the phase transition predicted by Patel's model is of second order, independent of 
the value of A^^- This originates from the fact that the 3D pure XY model has a second order 



phase transition although some subtleties with this result exist in the literature flS 



( see also |T^). The universality arguments of [|T4| on the other hand show that the QCD 
phase transition for A^^ = 3 is of first order, even in the case of pure gauge fields. Therefore, 
to resolve this apparent conflict, it is natural to include these vertices in these models and 
require them to be gluonic in nature. 

At first sight, Patel's results might seem a little unsettling as they rely on a mean field 
treatment as one is aware of the role of fluctuations in lower dimensions. However, as we 
will show below, this result is correct provided the results are interpreted in terms of large 
Nc- This is the issue we would like to address in this paper. 

We take the following route. The model that Patel gets on the lattice is the non-linear 
U(l) a model in three dimensions with some specific type of symmetry breaking. This 
model has a physical cutoff built into it, namely the thickness of the string. We evaluate 
the one-loop effective potential for this model in the continuum. We find that there are 
contributions which are non-analytic in the coupling constants. Being non-analytic they 
can't be absorbed by any finite choice of counter-terms and hence must be included in the 
initial theory. Therefore, we include the non-analytic terms in the free-energy expression 
and hence perform an "improved" mean-field theory calculation. 

The paper is arranged as below. In section 2, we discuss the behavior of these baryonic 
vertices with increasing Nc- Next in section 3, we discuss the fiux-string model for the 
simplest case of A^^^ = 2 where all hadrons can be represented simply as unbranched flux 
tubes and correspondingly the problem can be mapped into a non-backtracking random walk 
on a cubic lattice. We also review the model with A"c > 3 where the flux-tubes can branch 
and this warrants a different treatment from the = 2 case. In section 3, we give details of 
the mean fleld calculation done by Patel. In section 4, we discuss the naive continuum limit 
for the model of section 3 and discuss the phase structure using the continuum model. In 
section 5, we calculate the one- loop effective potential using the Coleman- Weinberg results 
and flnd the non-analytic contribution to the tree level effective potential. The non-analytic 
piece is then included and the mean fleld free energy is computed for the modifled system. 
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II. BARYONIC VERTICES AT LARGE Nc 



The string model can also serve as a model for large Nc QCD. In both models baryons are 
qualitatively differnt from mesons. In string theory this difference is coded in the iV-string 



vertex. Witten |20] has argued that the mass associated with the vertex increase with N. 
We will incorporate this behavior in our picture. The flux tube generated in three color 

2 

QCD has a finite cross sectional area. This area is expected to shrink as ^ in the limit as 
goes to infinity, where a is the width of the string and Qq is a constant. Since in the string 
model of a baryon there are N strings converging on the vertex the total cross sectional 
area of all the strings is finite. Thus there is finite region surrounding the vertex, of radius 
R ~ ao, where the strings overlap ( see Fig. p. 

It is this region of finite volume which provides the common background interaction 
for all the quarks ( with the strings) that is the essence of Witten's large N picture of 
the baryons. According to Witten ||20| this background interaction is responsible for a 
contribution of 0{N) to the baryon mass. Since, this junction of finite area is present for 
Nc , we incorporate this contribution as the effective mass v of the iV-string junction and 
thus have v ^ N as N oo. 

It is quite interesting that the above picture for baryons can be merged nicely with the 
Skyrmion picture for baryons. Since these flux tubes repel each other ( like their type-II 
superconductor counterpart), it is only natural that at large Nc, baryons would look like 
hedgehogs with quarks stuck at the ends of each "needle". We believe this picture is quite 
parallel to the case when one can interpolate between the bag or string picture for mesons 
121. 



III. FLUX TUBE MODELS FOR DECONFINEMENT 
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A. SU(2) flux-strings : Random Walk 



We review the flux-string picture of deconfining transition by describing tlie simplest 
case, when the color gauge group is SU{2). For SU{2) gauge theory, flux tubes represent 
both mesons (made of a quark and an anti-quark) and baryons ( made of two quarks). On 
the lattice, these flux strings can be represented by locally non-backtracking random walks 
||15|| . So, they are different from self-avoiding random walks, as they can cross each other 
after several steps. The Hamiltonian for a flux tube of length / is given by 

H = aal (3.1) 

and the grand partition function is 

Z(/3) = 5]e-(^"'") (3.2) 

For non-back tracking motion on a c?- dimensional hypercubic lattice, there are {2D — 1) 
different directions that the flux tube can move to and hence for a flux tube of length / in 
lattice units one has (2D — 1)' different configurations when the ends of the flux tube is 
allowed to be free. Hence the grand partition function is Q 

Z{p) =Y,{2D- l)'e-'^"''^ = e-'t^''^-''^^^'^-^)] . (3.3) 

The partition function diverges, which is a typical signature for phase transitions, at the 
temperature - 

ksT, = ^""^ , . (3.4) 
ln(2rf- 1) ^ ' 

In fact, this is a second order phase transition. This model has recently been discussed by 



Kiskis also [22 



one can include a chemical potential for quarks with finite masses, though we will ignore it 
hereafter. 
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B. Case for iV^ > 3 



The more general model proposed by Patel [15,16] can be described as follows. In this 
model one can include flux-string branchings which are present in models with N > 2, and 
hence this is more general than the model presented above. 

As before, we work on a D-dimensional hypercubic lattice. The fluxtubes are described 
via the the variables n^^^ living on the /x-th link on the lattice site i which are valued in the 
set {0, ±1}. Also one has quarks occupying the lattice sites described by the variable pi G 
{0, ±1} and the baryonic vertices described by the variable qi G {0, ±1}. The Hamiltonian 
describing the system is given by 

H = T.WaY^<, + ^P' + ^i] (3-5) 

i M 

where v is the energy associated with each A^-legged vertex representing the baryons. How- 
ever, as in the gauge theory, one also has the Gauss law constraint 

- ^i-f^,f^) -Pi + Nqi = = (3.6) 
which we need to implement at each lattice site. Hence the partition function is 

Z(aa,/3) = 5:e-^^n^a.o (3.7) 

i 

The partition function can be evaluated by rewriting the delta function as 

Kfl = ^ r de,e'^^'^ (3.8) 



This allows us to express our partition function as 

'^i,l-L Pi Qi 

Carrying out the summation over the variables ni^^,pi, qi we get 

de,. 



■ J ^TT m ■ n - n ■ 



ZiP)=I[J ^n(l + 2e-'''^"cos(e*+/'-^*))(l + 2e-^"cos^,)(l + 2e-^'^cos(iV^^,)) (3.10) 
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However, we are interested in the limit (3 oo and hence it is expedient to replace terms 
like (1 + 2e~^"^cos^j) by e^^°^^^ where h = 2e~^™. The partition function is then given by 

Z{P) = J pjc;6i.e-^SM,.=°"(^»+''"^«)+'^E»™se,+pX;,™s(Af0,) (3.11) 

i 

where 

J = 2e-^"'^ h = 2e-^'^ p = 26"^" (3.12) 



When the baryon vertices are very massive, i.e. p — 0, the expression (|3.11|) reduces to the 



partition function of the D dimensional XY model with an external magnetic field. 



The mean field free energy can be calculated ( as in Appendix 1) and one gets [|T^], 



f3F < -[ln/o(.) + DJi^-My + (/, _ + (3.13) 

Let us now deal with the case A'c = 3 in detail. For small values of z it is sufficient to 
study the critical parameters ( p.l3| ) using a series expansion in z for the free energy 

PFmj. = -y + -(1 - DJ) - -{p - 3h) - ^{1 - —J) + Oiz'). (3.14) 

Using this and assuming small values for p and h one can easily find the critical parameters 
16| . When the quarks are infinitely massive ( i.e. h = 0) one finds 

2 

l-DJo = ^ (3.15) 
3o 

where Jo is the critical parameter at the which the system shows phase transition. When 
p = , one can see Jq = and the system shows a second order phase transition. 

However, when quarks have some finite mass, h and the critical parameters are 
found to be 

2 3 

l-DJo = — Kr = —- (3.16) 

24 216 ^ ' 

When h > her, the system shows no phase transition. 

The above mean field treatment shows that the model exhibits a second order phase 
transition which is the case for = 2 and this result seems in agreement with the results 
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of [Q. However, as one knows mean field treatment is not reliable in c? > 4, one should 
treat this result rather cautiously. 

Before we attempt to discuss the stability of the phase structure against quantum fluc- 
tuations, it is interesting to discuss the phase transition in this type of model for arbitrary 

Nr > 3. 



In the second paper of Patel |16| it was claimed that the order of the phase transition 
depends on the value of p. When /i = 0, it was argued that for low values of p the phase 
transition is of second order and at some high values of p it is of first order. This can be 
checked explicitly by plotting the mean field free energy given by (|3.13|) ( see Fig. 
showing the critical cases). Though this is, in principle, true, we don't consider the result 
reliable. 

The second vacua appears at large values of z and for p near its maximum value of 2. 
This critical value of p seems to increase as Nc increases ( see Figures). On the contrary large 
Nc argument imply that p actually decreases as increases. Thus for large Nc, implying 
small p we expect only a second order phase transition. Given a physically reasonable value 
of p, (for Nc = 3), p < 1, Patel's model seems to imply a second order phase transition for 
Nc > 4, in agreement with [Q. This result also supports the recent arguments of |^ that 
a second order phase transition for Nc> 4: would provide insight into the weakly first order 
QCD deconfinement transition for Nc = 3. In our context, A^c = 3 is somewhat special. 



IV. RELIABILITY OF MEAN FIELD RESULTS 

To check validity of the mean field results, we have to compare one-loop contribution 
to the free energy with the mean field result. However, as it stands the Hamiltonian is 
awkward for such a perturbative calculation, since the interaction terms are nonpolynomial. 
The way to get around this is to notice that one can rewrite the Hamiltonian as the nonlinear 
0(2) ~ U{1) chiral model. Defining U = e*^ 
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Y: J{UiUl^ + h.c.) + h{U, + Ul) + p{U^ + u. 



(4.1) 



Notice that the quark masses appears hke the chiral symmetry breaking term as in the 
standard chiral Lagrangians but we also have the extra term due to the baryonic vertices . 
We also have to impose the constraint 



UU^ = 1 



(4.2) 



At long wavelengths we can replace the first term in equation ([4.1|) by the derivative so that 
the Hamiltonian reads 



(3H = J (fx^ [-JVU ■ VU^ + h{U + U^) + p{U^ + U^^) 



(4.3) 



where we have scaled the various coupling constants as follows. 

We will be interpreting the above model as an effective theory with a short distance 
cutoff ao which is the thickness of the string. This cutoff will be present as the fluxtube 
description breaks down at distances shorter than the string thickness. In terms of Oq one 
can see p = h = \ and J = —. We will restore this cutoff later. 

Even when p,h = this model is nontrivial due to its nonlinear structure and shows a first 
order phase transition which is understood as a vortex loop condensation |2^. This model 
also has an interesting large behavior as can be understood as follows. The quantum 
contribution to effective action for this model involves (A^ — 2)th order loop and hence has 
a factor of h^'"^. As A'c — > oo, this goes to zero and hence the theory is expected to become 
classical. 

Hereafter we will restrict ourselves to the case A^ = 3 only. Writing 



1 



(4.4) 



we get after eliminating i/j 

(5H = 



2 1 - 



:|V0|^- (/i-3p)0-4p0^ 



(4.5) 
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Note that in terms of variables 6 , (p = cosO. The potential term 



y(0)^_(/,_3p)0 



(4.6) 



admits a local minimum only if 3p > /i, which is given by 



3p — h 



= —V 



Lets expand the "Hamiltonian" ([4.5|) about this vacuum, 



(4.7) 



= —V + a 



(4.8) 



where a is the fluctuation about the vacuum v. Then one has 



r J 

(3H = / d^x\-Apv^ + — -Va ■ Va - Upva'^ - 4pa^] 

J 2(1 — 1>"') 



(4.9) 



Rescaling the field a = \/^^jX o^ig can rewrite the "Hamiltonian" for the fluctuation (after 



dropping the constant piece) as 



9X 



(4.10) 



where 



9p + h 



h 



J \ 3p 

1 ,9p+h.3 



VW J 

Note that when p = f the fluctuations become massless. 



(4.11) 



V. COLEMAN- WEINBERG EFFECTIVE POTENTIAL 

Though the nature of phase transition can be found by looking at the classical potential, 
the presence of fluctuations might change the structure in lower dimensions , i.e. mean field 
calculations based on Gaussian fixed points are not reliable. To find the deviations, it is 
customary to look at the one-loop effective potential, as calculated first by Coleman and 
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Weinberg Given a classical potential V{(j)) , the effective potential up to one-loop is 
given by 

For our case , ^(0) is as in ( [4 .61) and the inverse propagator G^^{k) is given by so that 
one-loop contribution to the the effective potential reads, 

The finite contribution from this integral can be evaluated by the dimensional regularization 
method, details of which can be found in the appendix. The answer is 

h 



12n 

where 



[W{<pp (5.3) 



W{<P) ^ -^0(1 - 02) (5.4) 

Note that the one-loop contribution is real when is negative. In terms of the original 
variables namely 6, this one-loop contribution reads 

_4v/6-(£)tcost^sin3^ (5.5) 

TC J 

This expression is non-analytic and hence can not be absorbed by a suitable counter. Hence, 
we follow the following strategy - namely, include this term in the outset and hence redo 
the mean field calculation. Basically, we are then working with an effective action which 
includes the non-analytic term. Restoring the cutoff one sees that the relevant term to add 
in the lattice theory would be 

-4y6-(^)tcost^sin=^^. (5.6) 

We would like comment on this quantum correction in the large Nc limit. It is known that 
in large limit the thickness of the strings goes to zero with ag ~ -kr- So, naively 
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the contribution to this term increases with Nc. However, if v ^ then one sees that this 
correction goes to zero a.s Nc ^ oo. Hence, the theory becomes classical, which is consistent 
with our naive expectations. 

In terms of the lattice theory, the inclusion of this term adds the following term in ( |A6|) 



-4V6-(^)t(cost^sin3 0) (5.7) 

71 JttQ 

However, this contribution is non-analytic in ^ = cos6. Hence, we have to take a branch-cut 
along the positive 2;-axis and take the positive root above the cut and the negative root 
below the cut. Then, one can readily evaluate the above integral, details of which we will 
relegate to the appendix, as 

, n ,3 TT ,5 „ ,5 7 , 9 ^ ,9 11 , , 

As a series expansion in z, this reads 

32^/6h, p ,3^^^^^ ^ ^^io. 2 nnnoo.3 nnioo.4 , 



TT^ Ja, 



)5 [0.022 + 0.013^ - 0.0013z^ - 0.0022^^ - .00188z^ + Oiz")] (5.9) 







this term will now contribute to the mean field free energy given by ( |3.14| ). Defining A = 
^^^^(-^)^/^, the total free energy is given by 

h 1 — 3 J V — 3h 
f3F = -0.022A - z{- + .013A) + z'^{ — - — + .0013A) - .0022A) (5.10) 

_ 3(1-4£)_.00188A) + .... 

64 

The behavior is this free energy function is similar to the mean field result and one can show 
that the phase transition do not change. 



CONCLUSIONS 

In this paper, we have attempted to incorporate quantum fluctuations in a flux-tube 
model by looking at a continuum limit of the model. The effect of the quantum fluctuations 
are interesting though they do not change the nature of the deconfining phase transition. 
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However, to show this one has to use the fact that the "Y- vertices" has to scale with Nc- 
Only then is the picture consistent with large N^. The large hmit of Patel's model implies 
that for Nc> 'i-i the phase transition is of second order. 
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APPENDIX A: MEAN FIELD CALCULATION 



Here we reproduce the Free energy expression for the model using the mean field cal- 
culation. Remember that our model is defined in three dimensions and hence mean field 
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calculations are not reliable. Then the partition function can be written as 

Zo= [ Y[dS{e)^e-^^^-^'-^'+^^''^^^-^' (Al) 

To evaluate the partition function, we now resort to mean-field where the spins are assumed 
to produce a mean magnetic field z so that the partition function is 

i 

= (^e-'™s(6»j+p-ei)+/icos0i+pcos{Ar0j)^^^^ (A2) 

where In{z) is the Bessel function with imaginary arguments of A^th order and (■ ■ ■)m.f. is 
the expectation value evaluated using the mean field weighing factor, i.e. 



{A) 



M.F. 



Now we use the fact for any convex function one has the following inequality 



e^) >e 



(A3) 



(A4) 



Hence the partition function is bounded below by the following quantity 

Z{(3) > [/g(^^)]^e('^cos(6»,+fj-(9i)+(/i-2)cos6»i+scos{Ar6»0> ^^5^ 

Recall that the free energy expression is given by F = lim^r^oo ~Jn^^^ hence 



(3F < (3Fm.f. = - 
Now notice that (sin i 



ln[/o(z)] + {Jj2cos{ei+^ - 9i) + {h-z) cosOi + scos{Nei)) 



and thus 



pF <- [ln[27r/o(z)] + 3(cos^i)2 + (/i - z){cos9i) + s{cosN9i) 

The various averages can be carried out easily by noting 

/o^ rf^e^'^^^^cos A^^ 



(A6) 



(A7) 



< cos A^^ > 



2"" „a cos 9 



/o " dO En=-oc In{a) COS uO COS A^^ 



Jq'^ dO T,n=-oo ^n{oi) COS hO 



which gives us the result 



< cos N9 >-- 



(A8) 



(A9) 



Using this is (|A8|) , we get the expression (|3.13|) . 
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APPENDIX B: THE ONE-LOOP CONTRIBUTION TO THE FREE ENERGY 



In this appendix we give the details of our various computations. From (|5.2| ), one can 
see that thew one-loop contribution can be written as 

VM = \j (0[ln(fc^ - 24^0(1 - 0^)) - In(^l^)]. (Bl) 

The second term in the integrand can be dropped in dimensional regularization. Defining 

W{(l)) = -2450(1 - 02)^ one can find 

dVM _ hw'{4>) f d^k 1 _ nw'{<p) 1 /— - 

-W ~ J (2vr)3 F + (0) - T6^^^-2) V ^^'^^ ^^2^ 



in 



iy(0)iy'(0) 



with W'{(f)) = Integrating both sides and dropping an irrelevant constant, we get 

VM = -^[Wm'/' = -^(£)^/^{0(l - 0)}^/^ (B3) 



Putting in = cos 6 one gets the answer ( |5.3|) . 

Next we calculate the mean-field contribution due to Vi(0). Now, 



< cos^/^ e sm-^ e > = — 5- 



lo e 

y — T ( / de cos=^/2+m ^ sin ^ - / de cos^/2+'" e sin . 



(B4) 



The integrals above can be performed by noting that Vcos^ is not an analytic function and 
we have to introduce a branch-cut as stated previously. We choose V cos e = 1 for ^ = and 



Vcose = — 1 for 6' = 27r. Then we can readily evaluate 

1 00 m 1 1 

< cos3/2^sin3^ > = — Vt E ^{ (B5) 

nIo{z) ^omrm + 5/2 m + 9/2^ ^ ' 

^ [^iFi(5/2, 7/2, z) - ^iFi(9/2, 11/2, ;z)] 



which leads us to (WM- 
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FIGURES 
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FIG. 2. Free Energy as a function of z for iVc = 4, %>= 1.48, J = .32 
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FIG. 3. Free Energy as a function of z for Nc = 5, p= 1.8, J = .32 
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FIG. 4. Free Energy as a function of z for Nc = Q, p = 2, J = 1/3 
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